Spin separation in digital ferromagnetic heterostructures 
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In a study of the ferromagnetic phase of a multilayer digital ferromagnetic semiconductor in the 
mean-field and effective-mass approximations, we find the exchange interaction to have the dominant 
energy scale of the problem, effectively controlling the spatial distribution of the carrier spins in 
the digital ferromagnetic heterostructures. In the ferromagnetic phase, the majority and minority 
carriers tend to be in different regions of the space (spin separation). Hence, the charge distribution 
of carriers also changes noticeably from the ferromagnetic to the paramagnetic phase. An example 
of a design to exploit these phenomena is given. 



PACS number(s): 75.70.Cn, 75.50.Pp, 75.10.-b 
The research in proactive use of the spins of the carriers 
to add a new dimension to electronics starts a new area 
known as spintronics. [Q The recent discoveries [^|J|] of 
ferromagnetism with high Curie temperatures in a num- 
ber of conventional semiconductors doped with magnetic 
impurities hold promise for the implementation of spin- 
tronics in semiconductors. 

Inhomogeneously doped semiconductors, such as the 
p-n junction, play a crucial role in conventional elec- 
tronic devices. Their properties depend on the distri- 
bution of the itinerant carriers, which is governed by 
the Coulomb interaction of the carriers with the impu- 
rities and with other carriers. In this paper we study 
the charge and spin distributions of the itinerant carriers 
in semiconductors delta-doped with magnetic impurities, 
such as GaMnAs. |^,|| Our theory is within the mean- 
field, effective-mass and virtual-crystal approximations. 
In addition to the electrostatic forces, the itinerant car- 
riers have an exchange interaction with the magnetic im- 
purities. Our calculation shows that the magnetization 
of the high Mn concentration in the delta layers in the 
ferromagnetic phase gives rise to a spin-dependent po- 
tential experienced by the carriers comparable in order 
of magnitude to the charge potential of the delta layer. 
The effect of the spin-dependent potential on the inho- 
mogeneous spin distribution of the carriers was noted by 
Loureriro da Silva et al. Q in multilayered GaMnAs with 
5% magnetic impurities. By contrast, the delta doping 
with a nominal concentration per atomic plane of 25% 
to 50% Mn atoms gives rise to a qualitative different 
phenomenon of spin separation. This creates the pos- 
sibility of the magnetic control of the itinerant carriers 
by manipulating the magnetization of the Mn ions. In 
particular, the spin-dependent potential may be used to 
influence the spin dependence of the distribution of itin- 
erant carriers. The majority spin carriers accumulate in 
the region of the Mn layers whereas minority carriers are 
repelled from the Mn region. In the paramagnetic phase, 
the spin potential averages to zero and the magnetic in- 
fluence disappears. This prediction of large changes in 
the distribution of carriers between the paramagnetic and 



the ferromagnetic phases may be a cause of the anoma- 
lous Hall effect in the ferromagnetic phase observed in 
the delta-doped systems. [Q To show the potential of the 
spin separation for device applications, we give an exam- 
ple of how heterostructures may be designed to engender 
spin separation. 

This paper is organized as follows. A brief review of the 
relevant theoretical formalism precedes the application to 
establish the general principle of magnetic control and to 
show how it works in two specific systems. The first is 
the experimental system of multiple digital layers of ref- 
erences [^1^ . We find the phenomenon of spin separation 
in the ferromagnetic phase. Specifically, for appropriate 
carrier densities and interlayer distances, the minority 
carriers are located mainly in the interlayer spacer, with 
a very small overlap with the magnetic atoms whereas 
the majority carriers are mainly located in the layers of 
magnetic atoms. Second, we design a double GaAs quan- 
tum well with AlGaAs barriers and a delta layer of Mn 
in the middle of one well and a delta layer of an accep- 
tor (Be) in the middle of the other. When the GaMnAs 
well is changed from the paramagnetic phase to the fer- 
romagnetic phase, a net transfer of charge and spin from 
the Be well to the Mn well is produced, because of the 
appearance of the dominant magnetic interaction. As a 
result, there is a potential drop across the structure as 
well as a spin polarization in the Be well. 

Our calculations are based on the mean-field model, 
|§-|12|,|| with two types of spins: localized d electrons 
with magnetic moments of 5/2 Bohr magnetons and itin- 
erant carriers (holes because of the Mn and Be doping). 
There is a Heisenberg spin exchange between the itiner- 
ant carrier and the Mn d electron, which is responsible 
for the ferromagnetism. The hole energy bands are based 
on the effective mass approximation with the kinetic en- 
ergy given by the bands of the host GaAs. The effective 
potential induced by the dopants is the only change for 
the carrier subbands, as in the virtual crystal approxima- 
tion. The Coulomb interaction between the holes is taken 
into account in the Hartrec approximation. The in-plane 
energy dispersion of the holes subbands is described by 
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a parabolic model which neglects the spin-orbit interac- 
tion. The omission of the spin-orbit effect, which has also 
been used in references PJ^-[lOt , greatly simplifies the nu- 
merical computation. On the other hand, it leads to an 
overestimate of the polarization of the Mn spins induced 
by the interaction with the carriers. We correct for this 
by using a smaller value of the coupling constant J of the 
Heisenberg exchange which reproduces the Curie temper- 
ature obtained in our previous calculation including the 
spin-orbit interaction. Because the charge and spin 
distributions are averaged properties of the system, we 
argue that the qualitative aspects of the effects discussed 
in this paper would remain if the spin-orbit interaction 
is included. 

Both model systems studied here are taken to be trans- 
lationally invariant in the xy-plane normal to the growth 
axis. Thus, the Mn concentration cm{z) varies only along 
the z axis. We assume that, in addition to the magnetic 
impurities which act as acceptors in GaMnAs, there is a 
distribution of donors, Cc{z), which partially compensate 
the Mn acceptors, so that the total density of holes is 
smaller than the density of Mn, as is observed. [|j For 
simplicity, we assume that the spatial distributions of 
the donor and acceptor impurities are the same except 
for a multiplicative constant. From the charge neutrality 
condition, P+Cc = Cm, where P, Cc and Cm are the av- 
erage densities of holes, of the compensating impurities, 
and of the Mn impurities, respectively. 

In the mean-field approximation, the effect of the mag- 
netic impurities on the itinerant carriers is given by a 
spin-dependent potential: 



V.^^-^cm{z){M{z)) 



(1) 



where a = ±1 denotes the spin directions, J is the ex- 
change coupling constant between the itinerant carrier 
and the Mn spin, and (M(z)) stands for the local aver- 
age magnetization of Mn. The Mn polarization is pro- 
duced by the molecular field created by the spin polar- 
ization of the itinerant carriers as well as the external 
magnetic field. At temperature T and zero external mag- 
netic field, the Mn magnetization is given by the usual 
Brillouin function: 



M[z) = S Bs 



J p+{z)-p^{z) 
knT 2 



(2) 



where Per (2) is the spin-dependent density of the itinerant 
carriers. 

The effective Schrodinger equation for the holes has a 
self-consistent potential Ve{z) + Va{z), where Ve{z) is the 
electrostatic potential seen by the hole given by, 

— 14(2) = — — [ce(z) - cm{z) + p+{z) + p^{z)] (3) 

The hole eigenenergy is the sum of the quantized energy 
En^a from the motion in the z direction and the plane 



wave energy e^:^^ = ^m'^ from the inplane motion. The 
hole wave function is the product of the bound state en- 
velope function in z times the plane wave normal to the 
z axis. The spin-dependent hole density is given by: 
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(4) 



where fj, is the chemical potential. 

It is illuminating to estimate the strength of the spin 
dependent potential (|l|) for a single magnetic digital 
layer. Throughout the paper we take the hole effective 
mass tensor to be TO|| = O.lf and = 0.37 times the 
free electron mass and model the distribution of the Mn 
in a digital layer as Gaussian: 



cm{z) 



C 



M 
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(5) 



In a digital layer of Gao.5Mno.5As the total concentration 
ofMnisCAf = 3.13xf0i'* cm^^ ^ 3.13 nm^^. Transmis- 
sion electron microscopy experiments Q reveal that the 
Mn is spread over 2 or 3 atomic planes. This corresponds 
roughly to the Gaussian half-width of A = 0.5 nm. The 
exchange integral of J = 150 mev-nm'^ used in our pre- 
vious work [p^ to calculate the Curie temperature is re- 
duced to J = 100 mev-nm'^ here in order to compensate 
the absence of spin-orbit in the present calculation. The 
maximum of the spin-dependent potential in Eq. for 
saturated magnetization is 441.5 meV in comparison with 
the charge well depth of about 177 meV. The peak to 
valley splitting is twice the maximum value. In a diluted 
magnetic semiconductor quantum well, the magnetic po- 
tential is one order of magnitude smaller than in a digital 
layer. Hence, the magnetic control is more practical in 
digital layers than in quantum wells. 

The magnetic potential is attractive for the majority 
carriers and repulsive for the minority carriers. As the 
temperature is increased, the magnetization decreases 
and so does the magnetic potential, vanishing above the 
Curie temperature. Accordingly, the itinerant carrier 
density profiles change significantly as the temperature 
varies. Both the potentials and the density profiles are 
shown in Fig. |l| for a system of iV = 4 digital layers, 
with an interlayer separation Z? = 40 ML (GaAs mono- 
layers), or 11.3 nm, and a density of holes per layer 
p = 1.3 • 10^'^ cm"2. Similar results are obtained for 
systems between 1 and 10 layers. In the left panel of 
Fig. |l] we present the potential profiles at T =40 K (thick 
dashed line), when the system is paramagnetic. In the 
same panel we show the total potential for majority (at- 
tractive) and minority (repulsive) carriers for spin de- 
pendent potentials at T =5 K when the system is fer- 
romagnetic. The large magnitude of the spin-dependent 
part of the potential is apparent in the figure. The effect 
on the density profile of the itinerant carriers is shown 
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in the right panel. In the ferromagnetic case (thin sohd 
hne) the carriers tend to pile up in the magnetic layers. 
In the paramagnetic case (long dashed line) the carrier 
distribution is more spread out, with more carriers in the 
spacer layers relative to the ferromagnetic case. 
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FIG. 1. Left panel: spin-dependent potentials for major- 
ity (solid line) and minority (thin dot-dashed line) at T=5K 
(the ferromagnetic phase) and total potential at T=40 K (the 
paramagnetic phase) (thick dashed line). Right panel: total 
carrier density profiles for T=5K (solid line, the ferromagnetic 
case) for T=40 K (long-dashed line, the paramagnetic case). 

Figure ^ shows the resultant spin-dependent density 
profiles for two temperatures. In the ferromagnetic case, 
the distributions of majority and minority carriers are 
distinct. The majority carriers are localized in the mag- 
netic layers whereas the minority carriers are almost 
totally expelled from it. This "spin separation" phe- 
nomenon is, of course, absent in the paramagnetic case. 
We have checked that as the interlayer distance decreases 
or the hole density increases the spin separation dimin- 
ishes. Spin fluctuations of Mn beyond the mean-field 
approximation are a major source of spin-flip scattering. 
Hence, we expect the spin-flip scattering to be reduced 
when the spin separation is larger. We have also found a 
correlation between the absence (or presence) of the ob- 
served anomalous Hall effect in transport experiments |^] 
and the presence (or absence) of the spin separation in 
our calculations of the same sample conflgurations. The 
anticorrelation may provide a clue to construct in future 
a theory of the anomalous Hall effect based on the ab- 
sence of spin separation. 

Results of figure (|^) show how the exchange interaction 
can overcome the electrostatic interactions in the case 
of a digital layer. We propose here a heterostructure 
in which this effect is made more apparent and can be 
measured. It is a double quantum well of GaAs with 



GaAsAl barriers, delta-doped with one layer of Mn in 
the middle of the left well and one layer of the acceptor 
Be in the middle of the right well. For illustration but 
not necessity, we choose the density profile of the Be to 
be identical to that of the Mn ions plus its compensating 
charges. We take the density of holes per digital layer in 



either well of 1.3 x lO^^cm ^. 
has a total density of 3.13 x 10^'^ cm^^ 



distribution as in Eq. 
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FIG. 2. Spin dependent hole density profiles for T=5K 
majority carriers (solid), minority carriers (long dashed) and 
for T=40 K (long-dashed). 

Above the Gurie temperature the system is totally 
symmetric (by the construction of the model) and so is 
the charge distribution (shown in the middle column of 
Fig. At low temperatures the majority carriers with 
spin antiparallel to the Mn are attracted by the mag- 
netic layer. As a result, there is a spin-dependent charge 
transfer from the non-magnetic to the magnetic well (the 
left column of Fig. ||), which generates an internal elec- 
tric field and a potential drop. In the rightmost panel we 
plot both the potential drop and the Mn average magne- 
tization as a function of the temperature. The potential 
drops follows the average magnetization and its maxi- 
mum is close to 20 meV at zero temperature. Larger 
values can be obtained for smaller carrier densities, but 
this decreases the Curie temperature |l^. The charge 
imbalance is 1.2 x lO^^cm^^ or 4.6% of the total hole 
density. The charge transfer effect takes also place for 
other values of the carrier density, as well as for struc- 
tures with a different density of Be and Mn. As a result 
of the rearrangment of the carriers in the heterostruc- 
ture, as the system goes from the paramagnetic to the 
ferromagnetic phase, large changes in the in plane and 
vertical transport properties can be expected as well. 
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FIG. 3. Left upper/lower panels: spin-dependent poten- 
tials/density profiles for majority holes (solid line) and minor- 
ity holes (thin dot-dashed line) at T=5K (the ferromagnetic 
phase). The middle panels: spin-dependent potentials and 
density distributions at T=40 K (the paramagnetic phase) - 
same notation as above. Right panel: potential drop (solid 
line) and magnetization (dashed line) versus temperature. 

Recent Monte Carlo simulations for bulk ferromagnetic 
semiconductors |l3| , in which the the Mn magnetic mo- 
ment is treated as a classical vector of fixed norm, confirm 
that the carriers are more localized around the magnetic 
atoms in the ordered phase than in the paramagnetic 
phase. This indicates that the attraction of the carriers 
to the magnetic atoms is a cooperative effect which is ac- 
counted for in our mean field calculation. Tight binding 
calculations for a superlattice of digital GaAs/GaAsMn 
give very similar results to those obtained using an 
envelope function approach [T^ ] , validating the use of the 
effective mass approximation in delta-doped systems. A 
density functional calculation for a superlattice of digital 
GaAs/MnAs predicts a half- metallic behavior ||l^. Our 
results should be compared with a similar calculation in- 
cluding compensating impurities, 50% of Mn and larger 
interlayer distances. 

In conclusion, we have studied the spatial distribu- 
tion of the carrier spins in digital ferromagnetic het- 
erostructures. Our main results are: (i) the distribu- 
tion of carriers in digital ferromagnetic heterostructures 
of GaAsMn is largely controlled by the exchange interac- 
tion, which overcomes the electrostatic interaction, (ii) 
Large changes in the spatial distribution of carriers take 
place when ferromagnetism is switched off (for instance, 
by increasing the temperature above the Curie tempera- 
ture), (iii) In some situations the minority carriers are to- 
tally separated from the Mn layers and the majority car- 
riers. This might have consequences in transport prop- 
erties and the anomalous Hall effect. 
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